A new FSA approach for in situ y-ray spectroscopy 



A. Caciolli a ' b '*, M. Baldoncini c , G. P. Bezzon d , C. Broggini a , G. P. Buso d , I. Callegari b , T. Colonna b , G. Fiorentini c ' d ' e , E. 
Guastaldi b , F. Mantovani ce , G. Massa b , R. Menegazzo a , L. Mou\ C. Rossi Alvarez a , M. Shyti c ' e , A. Zanon d , G. Xhixha c ' e ' f 

a lstituto Nazionale di Fisica Nucleare (1NFN), Padova Section, Via Marzolo 8 - 35131 Padova, Italy. 
b Center for GeoTechnologies, Via Vetri Vecchi, 34 - 52027 San Giovanni Valdarno, Arezzo, Italy. 
c Physics Department, University ofFerrara, Via Saragat, 1 - 44100 Ferrara, Italy. 
d Istituto Nazionale di Fisica Nucleare (INFN), Legnaro National Laboratory, Via dellUniversita, 2 - 35020 Legnaro, Padova, Italy. 
e Istituto Nazionale di Fisica Nucleare (INFN), Ferrara, Via Saragat, 1 - 44100 Ferrara, Italy, 
f Faculty of Forestry Science, Agricultural University of Tirana, Koder Kamez - 1029 Tirana, Albania. 



O 

> 
O 

r-- 



V. 

o 

• 1— I 
Oh 



Abstract 

An increasing demand of environmental radioactivity monitoring comes both from the scientific community and from the society. 
This requires accurate, reliable and fast response preferably from portable radiation detectors. Thanks to recent improvements in 
the technology, y- spectroscopy with sodium iodide scintillators has been proved to be an excellent tool for in-situ measurements for 
the identification and quantitative determination of y-ray emitting radioisotopes, reducing time and costs. Both for geological and 
civil purposes not only 40 K, 238 U, and 232 Th have to be measured, but there is also a growing interest to determine the abundances 
of anthropic elements, like 137 Cs and 131 I, which are used to monitor the effect of nuclear accidents or other human activities. 

The Full Spectrum Analysis (FSA) approach has been chosen to analyze the y-spectra. The Non Negative Least Square (NNLS) 
and the energy calibration adjustment have been implemented in this method for the first time in order to correct the intrinsic 
problem related with the^ 2 minimization which could lead to artifacts and non physical results in the analysis. 

A new calibration procedure has been developed for the FSA method by using in situ y- spectra instead of calibration pad spectra. 
Finally, the new method has been validated by acquiring y- spectra with a 10. 16 cm x 10. 16 cm sodium iodide detector in 80 different 
sites in the Ombrone basin, in Tuscany. The results from the FSA method have been compared with the laboratory measurements by 
using HPGe detectors on soil samples collected in the different sites, showing a satisfactory agreement between them. In particular, 
the 137 Cs isotopes has been implemented in the analysis since it has been found not negligible during the in-situ measurements. 

Keywords: environmental radioactivity, FSA, soil, calibration, sodium iodide 
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1. Introduction 

In situ y-ray spectrometry with sodium iodide scintillators is 
a well developed and consolidated method for radioactive sur- 
vey ( |Nuccetelli|[2008i |Chiozzi et al~l|2000[ |TyIer|[2007l ) with a 
wide range of applications from mineral exploration to envi- 
ronmental radiation monitoring ( [IAEA-TECDOC-1363 |2003 



providing quantitative information especially about abun dances 
of principal natural radioisotopes, 40 K, 238 U and 232 Th ( I AEA- 
ITECDOC-566 |1990| ). 

The experiences of nuclear power plant accidents and atmo- 
spheric nuclear weapon tests taught us that the fallout of man- 
made radioisotopes ( 137 Cs, 134 Cs and 131 I) can affect wide por- 
tions of the planet. In order to employ such powerful techniques 
in this context the sensitivity and the quickness have to be im- 
proved ( |Tyler|2004| ). 

The conventional approach for studying the specific activity 
concentration of the three principal natural radioisotopes is to 
monitor broad spectral windows during the analysis ( |Verdoya| 
et al. |[2Q09| pesbarats and Killeen||l990| ). Generally, these 
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windows are chosen around the photopeaks of 40 K (1460 keV), 
of 214 Bi (1765 keV), and of 208 T1 (2614 keV). Since the 238 U 
and the 232 Th are not y-rays emitters their concentrations are 
evaluated detecting the y-rays produced by 214 Bi and 208 T1 re- 
spectively. The assumption of secular equilibrium of the decay 
chains is required in order to use this approach. In addition to 
the abovementioned radionuclides, the three-windo ws method 
has been extended to the measurement of 137 Cs ( [Cresswell et 
|al.|2006t [Sanderson et al.[T989] ). 

The limit of this technique is that it becomes imprecise for 
short time acquisitions and it suffers the poor intrinsic energetic 
resolution of Nal(Tl) detector. In particular, the Compton con- 
tinuum around 137 Cs photopeak is generally very intense com- 
pared to the intensity of 662 keV y-line. This translates into 
long acquisition times. Moreover, the intrinsic energetic reso- 
lution of sodium iodide detectors prevents to resolve the triplet 
at 583 keV ( 208 T1), 609 keV ( 214 Bi), and 662 keV ( 137 Cs). This 
issue can be solved only by post processing the data. In any 
case the windows analysis method requires a prior knowledge 
of the presence on site of such radioisotope. As a matter of fact, 
this method is blind to unexpected radionuclides. 

Significant improvements in y-ray spectrum analysis have 
been obtained by implementing the full spectrum analysis 
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(FSA) method ([Hendriks et al.|2001| [Mnty|1992|[Gmllot|2001 



Gutierrez et al.||2002[ |Tyler||2Q07| ). Since the FSA uses the 
full extent of the spectrum, as opposed to the three windows 
method, it needs a much lesser statistic to reach the necessary 
accuracy. This, in turn, means a drastic reduction in acquisition 
times. 

In this paper the non-negative least square (NNLS) constraint 
and the energy calibration adjustment have been implemented 
into the FSA algorithm. This new approach to the FSA was ap- 
plied to a 10.16 cm x 10.16 cm portable Nal(Tl) detector dedi- 
cated to in situ y-ray spectrometry. 

Moreover, we experimented a new approach to the efficiency 
calibration procedure. Instead of using calibration pads, we 
used the spectra acquired in sites selected specifically for cali- 
bration purposes as described later in this article. 

A detailed characterization of each site (both for calibration 
and for measurement) was performed by sampling the soil in an 
area of about ten meters radius and by measuring, in the labora- 
tory, the radioisotopes concentrations by using a high resolution 
apparatus made of two lead shielded HPGe detectors (called 



MCA-Rad and described by |Xhixha et al.|2011| ). The Nal(Tl) 
detector was used within a larger project devoted to the inves- 
tigation of the soil characteristics in Ombrone basin (located 
in Tuscany Region) in morphologic, pedologic, and lithostrati- 
graphic way. About 80 different sites have been measured with 
two main finalities: the radioactivity characterization of soil 
for the study of lithological differences of the altered substra- 
tum and the validation of our method. The results of the new 
FSA analysis are furthermore compared to the characterization 
performed by measurements on sample in laboratory with the 
MCA-Rad apparatus and to the traditional three windows meth- 
ods. 



2. Calibration procedure and developments 

Portable instruments are usually calibrated by means of stan- 
dard spectra acquired at least using three concrete pads enriched 
in K, U and Th and a pad free of radioactivity representing the 
background (Hendriks et al. 2001; Engelbrecht and Shwaiger 
20081 |L0vborg et al.||1981 



These pads are usually concrete 
constructions of cylindrical shape, having finite dimensions of 
2-3 m in diameter and 0.3-0.5 m thick and for this reason it is 
needed a geometrical correction due to the differences from an 
infinite calibration sources. 

Hendriks et al. proposed as an alternative solution to build 
pads by using not concrete but silversand mixed with KC1 for 
the potassium pad, monazite for the thorium pad, and slags for 
the uranium pad. 

The design of an ideal pad with one radionuclide inside 
and with a perfect homogeneous distribution of the radioiso- 
tope in its volume is impossible, because impurities and non- 
homogeneities are always present. For example, a clear con- 
tamination of uranium in thorium pads has been often seen, as 
reported by Hendriks and co-workers. For this reason, compro- 
mises between accuracy and applicability of the method have 
to be weighted. In the case of in situ y-ray spectrometry accu- 



racies of the order of less than 15% are usually well accepted 
which legitimize the above assumptions. 

It is worth mentioning that the hypothesis of homogeneous 
distribution of the radionuclides in the pad mixture should be 
verified and that the cost of production, handling and storage of 
the pad is not negligible. 

Instead of building pads, an alternative calibration procedure 
is described in the present work. It is based on the selection of 
sites characterized by a prevalent concentration of one of the 
natural radionuclides. Even if it is almost impossible to select 
one site which contains only one of the nuclides, the selection 
will be oriented toward sites with unbalanced contents. Since 
most of the y-rays emanating from the earths surface originate 
in the top 50 cm ( |IAEA-TECDOC-1363 ||2003] ) depending on 
the rock/soil density the radioactivity characterization of each 
site is performed sampling at 10 cm depth. All calibration sites 
were selected using geological and geomorphological consid- 
erations and further validated trough laboratory measurements 
(a list of all the sites and the concentration of each isotope are 
reported in Table [TJ. The in situ measurements can be affected 
by the specificity of the place, like the soil non-homogeneity, 
the ground morphology, the non secular equilibrium in radioac- 
tive chains, the vertical distribution of 137 Cs, the presence of 
vegetation, moisture, etc. Thus, the calibration sites should be 
selected according precise prescriptions: 

• relatively uniform distribution of radionuclides in secular 
equilibrium with their products, 

• plane area without any steps and large enough to be ap- 
proximated as an infinite source (maximum 10 m radius), 

• undisturbed areas: assuring relatively constant 137 Cs verti- 
cal distribution 

• uniform and relatively homogeneous soil moisture content 
and vegetable coverage. 

For each site, a variable number of samples from 5 to 12 was 
collected in random positions within 10 m radius in order to 
check also the homogeneity of the site around the detector since 
more then 90% of y-rays detected by the sodium iodide are pro- 
duced by the 7 m radius and 0.5 m deep area around ( [Grasty 
et al. || 1979 ). The homogeneity of the each site is assured in- 
side the error reported in the in Table 1 Soil and rock sam- 
ples were dried, homogenized (less than 2 mm fine grain size) 
and sealed in measurement containers for at least four weeks 
in order to allow the 222 Rn growth up prior to be measured us- 
ing high-resolution y-ray spectrometry system MCA-Rad fol- 



lowing international standards of analysis ( |ASTM C 1402-04 
20091 |UNI 10797:|^999] ). Each sample is inserted between 



two HPGe detectors which are shielded with few cm copper 
and 10 cm lead in order to reduce the environmental radioac- 
tivity from the laboratory hall. Nitrogen gas is also fluxed into 
the lead shielding in order to remove residual radon gas. The 
MCA-Rad has been calibrated by using both certified radioac- 
tive sources (i.e. 152 Eu, 56 Co) and certified reference materials 
(RGK-1, RGU-1, and RGTh-1 produced by IAEA) containing 
known amount of natural radioisotopes of potassium, thorium, 
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Site 


K[%] 


eU [ppm] 


eTh [ppm] 


CA1 


53.9+0.7 


< 1.0 


6.0+0.5* 


CC2 


0.06+0.02 


0.7+0.3 


0.6+0.7 


GC1 


0.07+0.04 


0.27+0.08 


1.14+0.13 


GV1 


4.9+0.6 


7.5+1.1 


7+1 


PM2 


2.26+0.05 


2.27+0.18 


10.7+0.8 


RT1 


0.10+0.01 


6.8 +1.1 


1.74+0.16 


SMI 


1.54+0.14 


1.6 +0.3 


8.6+0.6 


SP2 


2.92+0.08 


7.5+0.4 


39+2 


ST2 


7.8+0.9 


36+5 


124+16 



a the Th is obtained by the procedure explained in section 3 



Table 1 : the average of the distribution of natural radioisotopes concentration. 
The errors correspond to one standard deviation. The conversion factors from 
Bq/kg are obtained by |IAEA-TECDOC-1363 |2003| 1 % = 3 1 3 Bq/kg for potas- 
sium, 1 ppm = 13.25 Bq/kg for uranium and 1 ppm = 4.06 Bq/kg for thorium. 



and uranium. These two different calibration procedures are in 
agreement within their systematics (about 5%). Only the CA1 
site (discussed in the next section in details) is not a outdoor 
site, but it is made by a pad of KC1 fertilizer. The concentra- 
tions reported in Table [T] are the average of the measurements 
performed in the laboratory on the collected rock samples. This 
way the heterogeneity of each site is properly implemented by 
the errors which are dominated by the spread in the results of 
the collected samples in each calibration site. The choice of a 
number of calibration sites greater than the number of analyzed 
elements, as in our approach, is a mandatory to avoid artifacts 
in the sensitive spectra. 
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Figure 1: spectra acquired with the sodium iodide in situ (black dashed line) 
compared with the fit obtained by the FSA+NNLS (see section [TT) algorithm 
(full line red). The 137 Cs contribution is shown alone (green dotted line) to 
underline the need to include this element in the analysis. 

One advantage of this calibration approach is the possibility 
to take into account radionuclides other than 40 K, 232 Th, and 
238 U (in general a minimum of k + 1 sites are needed, where k is 
the number of nuclides we want to be sensitive to). In our case 
a site with a prevalent concentration of 137 Cs has been included 
and used in the calibration (the CC2 site in Table[T]). As a matter 



of fact, after the Chernobyl accident, it is a mandatory to take 
into account the 137 Cs signal in the spectrum analysis. 

This is clearly visible in Figure [T] where the presence of ce- 
sium is definitively not negligible (the y-peak energy is 662 
keV). If necessary, other radioisotopes due to nuclear accidents, 
like 131 1, could be added in the future. 



3. The FSA Algorithm 

The full spectrum analysis method has been developed in dif- 

[20091 |Hendriks et"aT][200Tt 



ferent approaches (Maucec et al. 

|Minty|[llJ92l jCrossley and Reid||1982| [Smith et al.||1983| ) and 



was found to be a successful spectrum analysis tool. When ap- 
plied to scintillation detectors, it enhances their potentialities 
for radioisotope quantitative determination. As a matter of fact, 
the FSA encompasses almost the full energy spectrum, reduc- 
ing the required statistic of a single measurement and therefore 
its duration in time. 

The spectra acquired in situ are fitted by a linear combination 
of the fundamental spectra derived for each isotope from the 
calibration analysis. The events registered in each channel in 
the measured spectrum, N, can be expressed as: 



N(i) = Y d C k S k (i) + B(i) 



(1) 



k=\ 



where 

• N(i) are the counts in the channel i, 

• Cfc are the concentration of the element k, 

• Sk(i) are the associated counts to the fundamental spec- 
trum of the element k in the channel i, 

• B(i) are the counts in the channel i due to the intrinsic 
background. 

and the index k stays for 40 K, 232 Th, 238 U, and 137 Cs. It is be- 
come a conventional representation for in- situ measurements, 
for geological purposes, to express the concentrations of natural 
radioisotopes in their respective abundances, where 40 K is given 
in % weight while eU and eTh are given in ppm. The 137 Cs is 
expressed as the absorbed dose by the detector (nGy/h) due to 
the heterogeneous distribution property of anthropic radioiso- 



topes (described in section 3.2 ). The activity concentrations are 
deduced applying the least square algorithm to rectangular ma- 
trix and minimizing the^ 2 as in the following equation: 



X 2 = — = ~ E CkSk({) ~ *»] 2 /#(0, 

11 i=l k=\ 



(2) 



where the N(i) is considered poisson distributed and n is the 
number of channels in the spectrum used in the analysis (867 in 
our case). 

During the calibration of the system the fundamental spectra 
(the S matrix) are obtained by solving equation [T] with the ra- 
dionuclide concentrations (the Q coefficients) reported in Table 

m 



[S] = [C]- l x[N]. 



(3) 
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It has to be noted that FS A calibration method produces also the 
intrinsic background, B, which can be compared with a spec- 
trum acquired with the detector inside a thick lead shielding. 

Only the energy range from 300 keV up to 2900 keV is con- 
sidered in the analysis. Below 300 keV there is a strong pres- 
ence of the backscattering events which depends on the atomic 
number and density of the surrounding materials. Above 2900 
keV only the cosmic ray contribution is presenf] 

Once the first solution has been obtained, in order to im- 
prove thex 2 minimization, a trimming procedure is executed by 
changing the site concentrations (Q) in small steps around the 
measured intervals and repeating the matrix solution. This strat- 
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Figure 2: the spectra acquired in CA1 site in (red full line). The fit obtained 
using the concentrations measured with the MCA-Rad system is also reported 
(green dashed line). 



egy has been implemented in order to correct the unavoidable 
differences between the sample measured concentrations and 
the average site values. The need of this correction is evident 
observing the spectra acquired at the CA1 site where the sodium 
iodide detector was placed on a large platform made of KC1 fer- 
tilizer used as potassium calibration site. The KC1 was stored 
inside a building and in the Nal(Tl) spectra are clearly visible 
the peaks due to the thorium and uranium contained inside the 
building walls. In the samples analyzed by the MCA-Rad sys- 
tem, the contribution from these two elements was not present. 
Figure [2] shows the y-spectrum acquired in CA1 site compared 
with the fit spectrum obtained by assuming only the concentra- 
tions obtained by the measurements on samples in laboratory. 
Except for CA1, this procedure applied to the other sites pro- 
duce a deviation from the concentrations reported in Table [T] 
within lcr which confirms the assumption of acceptable homo- 
geneity declared in the previous section. 



lr The cosmic background has been found to be almost constant in our mea- 
surements and negligible with respect to the intrinsic background of the Nal(Tl) 
detector. However, its contribution is averaged and taken into account by the B 
sensitive spectrum determined by the algorithm. 



3.1. Improvements in FSA 

Thex 2 minimization without any further conditions, which is 
the base of the FSA method, can bring to sensitive spectra hav- 
ing energy regions of negative counts. Two evident examples 
of this problem are shown in Figure [3] and [4] 




1500 
Ey [keV] 



Figure 3: the sensitive spectra of 137 Cs, obtained using the standard FSA. The 
green line is placed to show the zero counts level. 

The presence of these non physical results introduces 
crosstalk effects in the analysis, leading to systematic errors. 



The NNLS (Non Negative Least Square) constraint (Lawson 
and Hanson] 1995; Desesquelles et al. 2009 1 Boutsidis and 



Drine a^|2009 ), which forces the counts on each bin to be zero 
or positive, has been for the first time implemented in the FSA 
algorithm in order to avoid this problem ( B aldon cmi] |20 1 0| ) . 
The NNLS algorithm is based on the Kuhn-Tucker condition 
for Problem of Least Square with Linear Inequality Constraint 
(LSI) ( [Lawson and Hanson 1995). In the case a problem has 
a solution with least square minimization this theorem states 
that there is also a solution of the problem with non negative 
solution. 

The sensitive spectra calculated with the new algorithm are 
shown in Figure [5] where it can be directly seen a more reliable 
sensitive spectra with the NNLS implementation. 

It has to be noted that the sensitive spectra resolution can be 
severely affected by gain mismatch between the different site 
spectra, which has to be always calibrate properly. Finally, 
residual correlations between isotopes are still present in the 
sensitive spectra shape, due to the presence of all the radionu- 
clides in most of the selected sites. 

3.2. l31 Cs 

The 137 Cs has been added to the analysis since it was found in 
not negligible quantities in soils. It was produced by the Cher- 
nobyl accident and randomly deposited by climatic events in 
Italy regions. The profile distribution of 137 Cs in soil systems 
shows a tenden cy of decrease toward deeper layers ( |Zhiyan-| 
ski et al. |2008 ). The 137 Cs is also heterogeneous in soils and 
this is confirmed by the measurements performed on collected 
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Figure 5: the sensitive spectra obtained through the FSA with NNLS constraint. 
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Figure 4: the sensitive spectra of 238 U, obtained using the standard FSA. The 
region where there are negative counts is emphasized in the box. 



samples which show a large distribution of 137 Cs values on the 
surface (see Table [2}. As a matter of fact, the Nal(Tl) detector, 
which records the y-rays produced by a wide surface area, is 
the best solution in order to average the cesium amount in the 
soil. In particular, it is possible to avoid over/under estimations 
due to point-like sampling. It has to be noted that the sensi- 
tive spectrum for 137 Cs is dominated by the CC2 site which is 
poor of natural radioisotopes and with an high value of 137 Cs. 
This way the cesium sensitive spectrum is free from correla- 
tions with the other radioisotopes. The dose absorbed by the 
detector, D, which takes into account the heterogeneity of the 
isotope in soils, is calculated summing the counts in each chan- 



site 


Bq/kg 


site 


Bq/kg 


CC2 


8 + 4 


GC1 


0.87 + 0.04 


RT1 


6 + 5 


SMI 


26 +37 

ZO -26 


ST2 


61+ ioo 

u -61 


SP2 


23 + 2 


GV1 


31 + 18 


PM2 


18 + 9 



Table 2: 137 Cs concentrations recorded on the collected samples in each cal- 
ibration site. The results of the samples are averaged and the errors are lcr 
standard deviation. 



nel of the sensitive spectrum S c s (i) weighted with the channel's 
energy E(i) and detector mass, m. This is multiplied by the Cq s 
coefficient derived by the FSA+NNLS algorithm as in the fol- 
lowing equation: 



D=-C Cs YSc s (i)E(t), 



(4) 



As described in Section |2| the sensitive spectra are calculated 
in a range of energy from 300 kev up to 2900 keV. The con- 
tribution to the dose due to energies below 300 keV has been 
included a posteriori by measuring the signal of a 137 Cs radioac- 
tive source and it brings a 5% systematic uncertainty. 

4. Spectra Analysis and Method Validation 

Recently an intense campaign of measurements was dedi- 
cated to the geological study of the soils in Tuscany, in an area 
investigated is situated between the regions of Siena and Gros- 
seto in the Ombrone's basin (see Figure [6]). This territory is 
characterized by a highly diversified geological structure and it 
has been studied by measuring different aspects and character- 
istics of the soil. The natural radioactivity of this area has been 
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also measured over 80 sites, where it has been investigated both 
by acquiring spectra in situ with the Nal(Tl) and by collect- 
ing samples for measurements with the MCA-Rad system. The 
sites selection fulfilled the same characteristics chosen for the 
calibration procedure, in particular for the topology of the site. 
For each site the 10.16 cm x 10.16 cm Nal(Tl) detector was 
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Figure 6: the position of the 80 sites used for the validation (red circles) in the 
area of Ombrone basin. 



used to perform 5 minutes ground measurements in- situ: the 
acquisition time was chosen in order to have enough statistics 
for a comparison with the three widow method too. Each spec- 
tra was analyzed using FSA method with and without NNLS 
constrain (based on the new calibration approach described in 
section [2} together with three windows method. In each site 
were collected one sample in the position where was grounded 
the detector and four samples were collected on the side bisec- 
tors of a 2 m side square centered at the grounded position of 
the detector. Each sample is treated in the same way as the one 
collected for the efficiency calibration of the system. The av- 
erages of the results obtained form measurements in laboratory 
for each site was used as reference to be compared with the 
outputs of the two methods. 

In Table [3] the correlation factor, Q, which minimize the rel- 
ative dispersion is obtained by using the following equation: 





K 


U 


Th 


x l 


IAEA 


1.12+0.07 


1.11+0.10 


1.00+0.09 




FSA 


0.99+0.06 


0.78+0.14 


0.86+0.07 


1.22+0.08 


NNLS 


1.00+0.06 


0.82+0.13 


0.92+0.07 


1.08+0.08 


NNLS opt. 


1.06+0.06 


0.87+0.12 


0.94+0.07 


1.06+0.05 



Table 3: the Q coefficients averaged for all the data samples. For the IAEA 
results the^ 2 is not shown due to the absence of a fit procedure. 



where, MCA and Nal are referred to the radioisotopes concen- 
trations calculated in laboratory and in situ. There is an agree- 
ment within the uncertainties between the FSA and the window 
method, but it seems that the window method usually overesti- 
mates the concentration while both FSA versions go in the op- 
posite direction. For the FSA method three different algorithms 
have been compared and the results are reported in Table [3j 
the standard FSA described in literature, the data obtained im- 
plementing the NNLS in the algorithm, and the FSA+NNLS 
with the optimization procedure introduced in the calibration 
method. This way, it is possible to understand the effect on the 
resulting isotope concentrations for each step by using the x 2 
as reference. All results reported in Table [3] agree with a factor 
Q = 1, which guaranties the reliability of the method for all 
elements. The correlation for the uranium element is affected 
by the atmospheric radon concentration at the time of the in situ 
measurement, although this discrepancy is within the uncertain- 
ties. 

The energy calibration adjustment has been included in the 
analysis of the measured spectra in order to correct for possible 
gain mismatching between the measured spectrum and the sen- 
sitive one. This problem is minimized implementing this post 
process correction as proven by the reduction of the^ 2 reported 
in Table [5] As an example, in Figure |7]the measured 137 Cs con- 
centrations, determined using both the standard FSA and in the 
FSA+NNLS methods, are reported. When the cesium abun- 
dances are low, the standard FSA algorithm introduces negative 
concentrations in order to have a better fit of the spectrum shape 
(corresponding to the 137 Cs line there is also the 609 keV peak 
due to the 214 Bi). 

Introducing the NNLS and the gain drift, the quality of fits is 
improved by 10% in average as shown in column 5 in Table [3] 

The correlation between the in situ measurements, analyzed 
by the FSA+NNLS, and the MCA-Rad measurements for all 



y ~ ^ MCA 2 

1=1 1 



(5) 



80 sites are reported in Figure [8j[9j and[10| The uncertainty re- 
ported on the MCA-Rad is the standard deviation of the average 
calculated over the five collected samples. 

The uncertainties on the Q factors are used as systematic un- 
certainties on the concentration measurements with the sodium 
iodide, since they contain both the contribution from the non 
homogeneity of soil ground and the systematics due to the anal- 
ysis algorithm. The uncertainties found in this way are: 5% for 
the potassium, 14% for the uranium and 7% for thorium, which 
are smaller than the ones requested for outdoor in situ studies, 
for civil and also geological purposes. 
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Figure 7: the measured Cs dose in the 80 sites determined with the new 
FSA+NNLS approach (red circles) and with the standard FSA (black triangles). 
The new algorithm avoids the negative counting and it reduces at the same time 
the uncertainties. 
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Figure 8: potassium concentration measurements done with the MCA-Rad sys- 
tem and with the FSA+NNLS analysis are plotted together. The average of 
five soil samples with relative uncertainty is plotted for the MCA-Rad analy- 
sis. No errors are associated with the Nal(Tl) data. The linear correlation line 
associated with the £1 value is shown (full line) with 1 <x error (dashed line). 



5. Summary and Conclusions 

The in situ y-ray spectroscopy with sodium iodide scintilla- 
tors is one of the most powerful techniques for the measure- 
ment of radionuclides concentration. The main advantages are 
its modest time consumption, its portability, and the reasonable 



cost of the detector with respect to competing systems (Sai'dou 
|et al.|2008] ). These characteristics make the use of sodium io- 
dide scintillator a prime candidate to meet the recent increase 
in the demand for in situ y-ray spectroscopy. 

The FSA is a powerful tool for y-spectra analysis thanks to 
its reduction in required statistic and its increase in analyzable 
radionuclides. It is now also applied to airborne and in situ y- 
ray spectrometry. 



Figure 9: uranium concentration measurements done with the MCA-Rad sys- 
tem and with the FSA+NNLS analysis are plotted together. The average of 
five soil samples with relative uncertainty is plotted for the MCA-Rad analy- 
sis. No errors are associated with the Nal(Tl) data. The linear correlation line 
associated with the Q, value is shown (full line) with 1 cr error (dashed line). 
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Figure 10: thorium concentration measurements done with the MCA-Rad sys- 
tem and with the FSA+NNLS analysis are plotted together. The average of 
five soil samples with relative uncertainty is plotted for the MCA-Rad analy- 
sis. No errors are associated with the Nal(Tl) data. The linear correlation line 
associated with the Q, value is shown (full line) with 1 cr error (dashed line). 



Nevertheless, the FSA shows limits due to the intrinsic x 2 
minimization which brings to non physical results, both in the 
fundamental spectra construction and in the determination of 
the element concentrations. These problems have been solved 
in the present work by introducing the Non Negative Least 
Square constraint into the FSA algorithm. More reliable fun- 
damental spectra have been obtained and the uncertainties in 
the fit procedure have been reduced. 

At the same time, a new calibration procedure has been in- 
vestigated to reduce the difficulties connected with pads cali- 
bration. This way, 137 Cs isotope has been also introduced in the 
analysis taking into account the anthropic effect on the environ- 
ment which is clearly not negligible. 
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The new algorithm has been validated with measurements 
performed in 80 sites in Tuscany by using a 10.16 cm x 10.16 
cm Nal(Tl) detector. In particular the concentrations have been 
obtained with a 5% error on 40 K, 7% error on 232 Th, and 15% 
error on 238 U concentration. 

The FSA+NNLS method has been used also to analyze air- 
borne measurements where the need of advantage of working 
with reduced statistics allows to reduced equipment weight and 
flying costs. This analysis approach has been also found to 
be the optimal solution in order to characterize large amounts 
of different sites avoiding the necessary of sampling measure- 
ments reducing the analysis time and avoiding the treatment of 
samples. These results will be the subject of future publications. 
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